Revisiting two holes in a locally antiferromagnetic background: 
the role of retardation and Coulomb repulsion effects 
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The problem of two holes in the presence of strong antiferromagnetic fluctuations is revisited using 
computational techniques. Two-dimensional clusters and 2-leg ladders are studied with the Lanczos 
and Truncated Lanczos algorithms. Lattices with up to 2 x 16 and \/32 x \/32 sites are studied. 
The motivation of the paper is the recently discussed spatial distribution of holes in ladders where 
the maximum probability for the hole-hole distance is obtained at d = \/2 in units of the lattice 
spacing, a counter-intuitive result considering that the overall symmetry of the two-hole bound 
state is d^2_y2. Here this effect is shown to appear in small ladder clusters that can be addressed 
exactly, and also in planes. The probability distribution of hole distances d was found to be broad 
with several distances contributing appreciably to the wave function. The existence of holes in 
the same sublattice is argued to be a consequence of non-negligible retardation effects in the t — J 
model. Effective models with instantaneous interactions nevertheless capture the essence of the hole 
pairing process in the presence of short-range antiferromagnetic fluctuations (specially regarding the 
symmetry properties of the condensate), similarly as the (non-retarded) BCS model contains the 
basic features of the more complicated electron-phonon problem in low temperature superconductors. 
The existence of strong spin singlets in the region where the two hole bound state is located is here 
confirmed, and a simple explanation for its origin in the case of planes is proposed using the Neel 
state as a background, complementing previous explanations based on a spin liquid undoped state. 
It is predicted that these strong singlets should appear regardless of the long distance properties of 
the spin system under consideration, as long as the bound state is d^2_y2. In particular, it is shown 
that they are present in an Ising spin background. The time retardation in the family of t— J models 
leads naturally to low-energy hole states with nonzero momentum and spin one, providing a possible 
explanation for apparent SO(5)-symmetric features observed recently in this context. Finally, the 
influence of a short-range Coulombic repulsion is analyzed. Rough estimations suggest that at a 
distance of one lattice spacing this repulsion is larger than the exchange J. The hole distribution 
in the d^2_y2 bound state is reanalyzed in the presence of such repulsion. Very short hole-hole 
distances lose their relevance in the presence of a realistic hole-hole interaction. 

PACS numbers: 74.20.-z, 74.20. Mn, 74.25. Dw 



I. INTRODUCTION 

Among the most appealing scenarios to explain the 
behavior of the high temperature superconducting com- 
pounds are those based on antiferromagnetic fluctuations 
as the mechanism for hole pairing. These ideas have been 
formulated over the years in a variety of contexts. Early 
diagrammatic work for heavy fermions in the presence 
of antiferromagnetic correlations suggested that magnon 
interchange can lead to d-wave pairing The nearly 
antiferromagnetic Fermi liquid proposal arrives to the 
same conclusion also in a diagrammatic context where 
low-energy spin excitations are assumed to coexists with 
holes [|| . Since in real doped cuprates the antiferromag- 
netic correlation length S^af is only of a few lattice spac- 
ings other approaches have emphasized the relevance of 
the short distance physics in the problem Indepen- 
dently of these approaches which are based on the resum- 



mation of particular sets of Feynman diagrams, a vast 
amount of computational work has been devoted to the 
interaction of holes with spin excitations jj]. Although 
these techniques have as a natural limitation the finite 
size of the clusters that can be studied, they have the 
important advantage that inside these clusters the infor- 
mation gathered is basically exact. Since ^af is small in 
the cuprates, and, in addition, experimental estimations 
of the coherence length also suggest small size Cooper 
pairs 1^, computational approaches are very important 
to clarify the physics of carriers in a fluctuating antifer- 
romagnetic background. Following these ideas the exis- 
tence of pairing of two holes in the 6,^2 _y2 channel using 
the t — J and Hubbard models was observed since the 
early theoretical studies of the cuprates |4||^ . For realis- 
tic values of J /t the size of the pair has been known to be 
small, of only a couple of lattice spacings at most . In 
addition, recent work in the limit where the attraction 
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between holes is assumed much larger than the band- 
width provided a real-space simple qualitative picture 
for the formation of dr^2_y2 bound states It turns 
out that not only a suitable effective potential between 
holes is needed to obtain pairs in the d-channel (repulsive 
on site, and attractive between nearest-neighbors) but in 
addition the form of the dispersion of holes renormalized 
by spin fluctuations plays an essential role in stabilizing 
d-wave pairing instead of extended s-wave pairing Q . 

The small tight pairs found numerically suggest that 
the cuprates are in an intermediate regime between the 
BCS limit, where the pair size is much larger than the 
mean distance between carriers, and the Bose condensa- 
tion regime defined in the opposite limit. Recent experi- 
mental work reporting the presence of pseudogap features 
in the normal state of the underdoped cuprates |^ have 
revived the possibility that hole pairs may exist above 
the superconducting critical temperature, at least as fi- 
nite lifetime fluctuations |jl^ . These results provide extra 
support to the idea that a "real-space" approach to hole 
pair formation is more enlightening than a momentum- 
space approach. In the context of holes in antiferromag- 
nets it is natural to propose an effective model for holes 
where these particles move within the same sublattice (as 
suggested by a variety of studies of the one hole prob- 
lem in an antiferromagnet), interacting through an effec- 
tive potential induced by antiferromagnetism |11[ . In the 
limit where this potential is assumed to be short-ranged 
the presence of d-wave bound states has been deduced 
analytically through the solution of the two-body prob- 
lem, and the presence of superconducting correlations has 
been found using the BCS mean-field approach with RPA 
fluctuations supplemented by numerical studies |^ . More 
recent work has extended these ideas to include the in- 
terchange of magnons, leading to a longer range effective 
potential between holes reaching similar conclusions re- 
garding the d-wave pair formation in the problem 12|. 

Recently, a complementary approach to the study of 
holes in the cuprates has been considered. It is based 
on the analysis of even-leg ladder systems since they 
have a robust spin gap, a feature that is also present in 
the normal state of underdoped cuprates. The 2-leg lad- 
der has an antiferromagnetic correlation length of about 
3 lattice spacings iQ , similar to the correlation observed 
in the doped cuprates. Then, ladders have coexisting spin 
gap features and short-range antiferromagnetism. In this 
context it has been observed that a couple of holes intro- 
duced in the 2-leg ladder form a bound state for realistic 
values of J/t with characteristics of a d-wave pair [ p^Jl6| 
(although strictly speaking a sharp distinction between 
d- and s-wave pairs does not exist in ladders) . While the 
main reason for the hole attraction is the minimization 
of broken rung singlets [T7| , the appearance of a d-wave 
character in the pairing is likely caused by the presence 
of short-range antiferromagnetic fluctuations, establish- 
ing an interesting analogy between planes and ladders. 

Recent work in this context using the Density Matrix 
Renormalization Group (DMRG) |18[ approach has con- 



firmed the presence of two hole bound states in the d- 
channel in the case of 2-leg ladders jio). However, in this 
study it was remarked that holes in the bound state spend 
part of the time located along the diagonals of the ele- 
mentary plaquettes in the problem, a feature somewhat 
counter-intuitive for d-wave pairs since a two-body prob- 



lem with drr 



character cannot have particles in such 



a configuration. In addition, the presence of strong spin 
singlets in the vicinity of holes, resembling the Resonat- 
ing Valence Bond (RVB) states characteristics of spin 
liquids were noticed. In ladders an explanation for 
these results was proposed in Ref. : the singlets along 
plaquette diagonals have frustrating character and pair- 
ing occurs to share frustration. As a consequence of this 
effect domain walls are induced in the problem. A sim- 
ilar numerical result in the context of two-dimensional 
planes was earlier reported using Exact Diagonalization 
(ED) techniques on small clusters | |2l|] . 

The above mentioned recent numerical studies moti- 
vated in part the present paper. Its purpose is to fur- 
ther analyze the two hole problem in planes and 2-leg 
ladders using Exact Diagonalization [Q and Truncated 
Lanczos | |2^ ] algorithms, as well as to propose intuitive 
arguments to explain the results observed in the present 
and previous investigations. Thephysics obtained with 
the DMRG method on ladders is found to be con- 
tained also in small clusters that can be handled exactly. 
An intuitive picture both in real and momentum space is 
provided to justify how holes can exist in the same sub- 
lattice in a da;2_y2-wave bound state. The importance of 
retardation effects for this feature is remarked, and thus 
it is unavoidable to conclude that time-dependent hole- 
hole potentials arc needed to quantitatively account for 
the physics of the t — J model. However, effective mod- 
els with instantaneous interactions are expected to cap- 
ture most of the relevant physics, specially regarding the 
symmetry of the condensate, similarly as the BCS model 
does for the electron-phonon problem. The existence of 
strong spin singlets in the immediate vicinity of tightly 
bounded holes on planes is here argued to be caused, at 
least in part, by the d-character of the bound state and 
its existence is independent of the special features of the 
spin background. In particular they exist in the case of 
an Ising background, which does not have quantum fluc- 
tuations. For the case of ladders, since spin singlets are 
already formed in the undoped system our proposed ex- 
planation only addresses the unexpected strong strength 
of some particular spin singlets, and it is thus comple- 
mentary to previous explanations [l9| . It is also reported 
here that the special configuration where holes are lo- 
cated at short distances, such as V2 lattice spacings, is 
just one of several equally important hole distances in 
the two hole problem, and in addition it becomes unsta- 
ble after the introduction of a realistic nearest-neighbor 
(NN) Coulombic repulsion. 
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II. HOLES IN SAME SUBLATTICE IN BOUND 
STATES TRANSFORMING AS - Y'^ 

A. Discussion in Real Space 

As explained in the Introduction, previous stud- 



ies 



have emphasized that the configuration where 
the two holes are located in the same sublattice across the 
diagonals of a plaquette has a substantial weight in the 
two-hole ground state. At first sight this result seems 
counter-intuitive since one is used to the naive notion 
that in the d,j.2_y2 subspace the wave function of two par- 
ticles has a node along the lattice diagonals. However, 
remember that this picture is based on the idealization 
of the two-hole problem in a spin background (which is a 
A''— 2 body problem, where N is the number of sites of the 
cluster under consideration) as a system of only two par- 
ticles in an otherwise empty lattice interacting through 
a static effective potential. This potential is staggered in 
real space (peaked at (tt, tt) in momentum space) , favor- 
ing the location of holes in opposite sublattices. A wide 
variety of calculations including diagrammatic techniques 
supplemented by Quantum Monte Carlo studies |^ , i — J 
model estimations , and other approaches all agree 
in this respect. Using such a hole- hole potential there is 
no doubt that the two particles must be located in oppo- 
site sublattices for the d^2_y2 symmetry to be realized. 
Even within the context of the full t — J model, if it were 
possible to "integrate out" the spin degrees of freedom, 
the resulting potential between the two holes at zero fre- 
quency must again induce the same feature. 

For the case of planes, the key ingredient to understand 
the apparent discrepancy between the numerical and an- 
alytical calculations detailed above are the retardation 
effects caused by the finite velocity propagation of the 
spin excitations in lightly doped antiferromagnets. As a 
simple example consider a couple of holes located next to 
each other in a perfect Neel spin background, as shown 
in Figs.la-b. To construct a d.j.2_y2 wave function in this 
context, it is necessary to combine Fig. la and lb with 
a negative relative sign. This provides the simplest ex- 
ample in the two hole subspace of a contribution to the 
dx2_y2 bound state. If the Neel state is considered as 
the reference "vacuum", then each hole at site i can be 
formally assigned an effective spin according to the def- 
inition Sl^^ ^ Si — Sf^, where Si is zero for a hole 
and Sf^ is the spin of the reference state at the same 
site. With this definition holes at distance of one lattice 
spacing have opposite effective spins. However, there are 
other contributions to the many-body ground state of 
two- holes in the d^^-y^ subspace caused by the presence 
of a nontrivial spin background. For example, Fig.lc il- 
lustrates what occurs if a hole is moved by one lattice 
spacing from its position in Fig. la. In this situation the 
spin that is moved in the opposite direction of the moving 
hole is now located in the wrong sublattice with respect 
to the background. Once again if the Neel state is de- 



fined as the reference vacuum, such an arrangement of 
holes and spins is represented as in Fig. Id showing that 
it corresponds now to a three — body problem. The holes 
being now in the same sublattice have the same spin pro- 
jection according to the definition of the spin S^^^ given 
above, and the overall spin conservation is achieved by 
realizing that the actual spin involved in the hole hop- 
ping process, which now lives in the wrong sublattice, 
plays the role of a spin one excitation with respect to the 
Neel background. Thus, the hole hopping in the lightly 
doped t — J model can be rephrased as the destruction of 
a hole with spin a, and the creation of a hole with spin 
—a and a bosonic excitation of spin 2cr (Note that the 
pair operators introduced in Ref. |^ also contains the 
notion that spin excitations are needed to have d-wave 
symmetry in this context). Formally this corresponds to 
a hopping process regulated by the Hamiltonian 



hopp ' 



E 



(4+u,-^Ci,^a[2^ + h.c.) 



(1) 



where the c-operators represent the (fermionic) holes, 
and the a-operators the (bosonic) spin one excitations. 
The rest of the notation is standard. Once the three- 
body character of the problem is understood, it is easy 
to find combinations of individual states that transform 
globally as — y^- For instance, Fig.le contains four 
states with holes located along the diagonal of a plaque- 
tte coexisting with a spin excitation as described before. 
If the four states have equal absolute weight but the rela- 
tive signs indicated in Fig.l.e, then it can be shown that 
a dr^2_y2 state is formed. This reasoning has been tested 
by solving a 3 x 3 cluster exactly fixing one hole at the 
center and allowing for a second one to orbit around. 
Simply considering the dominant configurations in the 
ground state, the three-body characteristics of the mo- 
tion of two holes around each other become clear. In 
addition, certainly the discussion based on Figs.la-e can 
be extended to distances larger than those involved in a 
single plaquette. Actually as J/t is reduced in the two 
hole problem, longer and longer distances between the 
holes will have larger weights in the wave function since 
J/t regulates the strength of the attraction. 

Based on this reasoning the computational results pl| ] 
reporting that in planes the configuration with holes lo- 
cated in the same sublattice has the largest weight in the 
two hole ground state can be rephrased as an indication 
that the spin excitations are not "instantaneous" from 
the point of view of the holes, but instead they have a fi- 
nite lifetime. In other words, the process of holes moving 
within a given sublattice is not a rapid tunneling pro- 
cess. "Retardation" effects apparently are important in 
models where antiferromagnetic correlations at short dis- 
tances are robust. This is not too surprising considering 
that a typical spin wave velocity is regulated by the ex- 
change J while a variety of numerical results have shown 
that the quasiparticle bandwidths in the t— J and related 
models for cuprates are also of order J, at least in the 
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FIG. 1. (a-b) Pictorial representation of a couple of holes 
in a spin Neel state; (c) The movement of a hole in (a) intro- 
duces a spin incorrectly aligned with respect to the staggered 
order in the background; (d) Taking the Neel state as refer- 
ence (c) corresponds to a three body problem, namely two 
holes and a spin excitation; (e) The three bodies of (d) can 
be arranged in a d^2_y2 state as indicated. 

Note that the effect presented here is "local" in the 
sense that it involves just a few sites of the lattice. 
Thus, naming the spin excitation of Fig.lc as a spin- 
wave, which is usually associated to a delocahzed object, 
is somewhat misleading since a localized excitation can 
actually be decomposed into planes waves carrying a sim- 
ilar weight for all momenta. Since the density of states 
of magnons peaks at k = (7r/2, 7r/2), these states may 
matter more for local processes than the low-energy ex- 
tended excitations of momentum (7r,7r). This is similar 
to what occurs in electron-phonon problems where the 
phonons at the Debye frequency are more relevant than 
the low-energy acoustic modes. Then, the spin-wave ve- 
locity (which is defined in the vicinity of (tt, tt) for spin 
problems) is likely not the most important quantity to 



judge effects of retardation in the formation of tight hole 
pairs. 

B. Strong spin singlets along plaquette diagonals: an 
intuitive explanation 

The DMRG studies of Ref. ^ for the case of two 
holes located at distance a/2 also reported the existence 
of a strong spin-singlet along the diagonal the opposite 
to the one where the holes are located (for other recent 
studies see Ref. [^). In the context of ladders their 
existence was argued to be caused by the presence of 
frustrating singlet effects [|l9|. In addition, singlets are 
already formed in the ground state of 2-leg ladders and, 
thus, it is natural that they contribute to the movement 
of holes. However, the strong strength of these "diago- 
nal" singlets is unusual, and even more strange is that 
they also appear in two-dimensional clusters as shown 
numerically below. Then, the authors believe that there 
must be some other ingredient contributing to the exis- 
tence of these strong diagonal singlets. The discussion of 
the previous subsection allow us to propose an explana- 
tion based on holes in an antiferromagnetic background 
that complements previous discussions [^9| based on spin 
liquid states. Consider once again Figs.la-b remember- 
ing that these two states enter in the wave function of a 
dj.2_y2 state with opposite signs. Moving the right hole 
in Fig. la up one lattice spacing, Fig.lc is obtained where 
the spins along the diagonal opposite to the holes have 
opposite z-projections. This is once again reproduced in 
Fig. 2a. Now consider Fig. lb and move the upper hole 
to the right one lattice spacing: in this case (Fig. 2b) the 
spins along the diagonal opposite to the holes are again 
antiparallel, but with projections spin inverted with re- 
spect to those found in Fig. 2a. Next, combine Fig. 2a 
and 2b taking into account the fact that Figs.la-b have 
a weight of opposite sign in the ground state, and that 
the spins not belonging to the plaquette being analyzed 
have not changed. Then, it is clear that a spin singlet 
along the diagonal opposite to the holes is formed, as il- 
lustrated in Fig. 2c. Note that this reasoning should work 
specially well in an Ising background (i.e. using a Neel 
state without quantum fluctuations), and in Sec. Ill nu- 
merical results confirm this prediction using the t — 
model ||2^. Note also that for the case of an extended s- 
wave, the same line of arguments would lead to a triplet 
along the diagonal, a prediction that will be tested nu- 
merically below. 

A similar reasoning helps in understanding why sin- 
glets are formed next to hole pairs even if they are not 
along diagonals. Consider Fig. lb and move the upper 
hole first to the right and then down one lattice spac- 
ing. The resulting spin configuration is shown in Fig. 2d. 
Combining this state with Fig. la and, again, remember- 
ing that Figs, la and lb had a weight of opposite sign, now 
Fig.2e is obtained. Then, the presence of singlets next to 
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holes is also natural when d.j.2_y2 states are considered, 
and they are there regardless of the overall ground state 
properties of the spin system. 
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FIG. 2. Spin and hole configurations allowing us to show 
that strong spin singlets are expected in the vicinity of holes 
orbiting in a d^2_y2 state when a Neel background is consid- 
ered. For details see the text. 



C. Discussion in momentum space 

The previous discussion carried out in real-space in the 

limit where holes are heavy can be reinterpreted using 
a diagrammatic approach in momentum-space consider- 
ing fermions (holes) interacting through the exchange of 
bosons (spin excitations). Fig. 3a contains a " randomly" 
selected diagram contributing to the bound state of two 
holes in a spin background. The important point to re- 
mark arises if a snapshot of the system at an arbitrary 
time (such as the one indicated by the dashed line) is 
taken: in this case it is clear that together with the pair 
of fermions there may be a nonzero number of spin excita- 
tions. In spite of this increase in the number of particles 
at intermediate times, all must be combined in such a 
way that the original spin and symmetry under spatial 
rotations of the fermionic pair at time t = — oo is con- 
served. In the case of relevance for the two-hole bound 
state the total momentum and spin must remain zero, 
and the d^^-y^ symmetry under rotations must be pre- 



served. Consider the intermediate time denoted by the 
dashed line in Fig. 3a: in this situation one vertex interac- 
tion has taken place at previous times, changing the spin 
of the fermion involved, and now both fermions have the 
same spin projection with the boson carrying the com- 
pensating spin, as in the real-space picture of Fig. Id. In 
addition, if the boson has momentum q, the holes orig- 
inally in — k and k, now switch to — k and k — q. At 
this intermediate time and if only the holes are consid- 
ered one may arrive to the counter-intuitive conclusion 
that they have switched to a state of spin 1 and momen- 
tum q, that in addition is not in the d,j.2_y2 subspace (in 
principle they can be in any irreducible representations 
of the Civ symmetry group of the square lattice). Such 
a conclusion would be incorrect since the overall quan- 
tum numbers are given by the combination of fermions 
and bosons. Concentrating on the behavior of only the 
fermions in the problem may lead to apparent paradoxes, 
as described before. 





(a) (b) 

FIG. 3. (a) Momentum space representation of the inter- 
action of two holes (solid lines) mediated by spin-waves (wavy 
lines). Spin and momenta arc indicated. The dashed lino rep- 
resents an intermediate time relevant for the discussion in the 
text. Time runs upwards in the graph; (b) The same Feynman 
diagram as in (a) but now assuming that the interchange of 
spin- waves is instantaneous (no retardation). This collapses 
the spin- wave propagators into a point. Momenta and spin 
are indicated. 

If the bosonic mediated interaction is considered in- 
stantaneous, i.e. when a w = potential is used to de- 
scribe the hole-hole interaction, then a situation where 
at all times the holes carry the quantum numbers they 
had sX t = —00 is recovered (see Fig. 3b). In this case, 
holes miist necessarily reside; on opposite sublattices for 
its bound state to transform as x^—y"^. Then, it is the fact 
that the w ^ Q effects (retardation) cannot be neglected 
in the t — J model that allows for hole configurations be- 
longing to the same sublattice. The diagrams in Figs. 3a 
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and 3b are similar to those arising in the context of low 
temperature superconductors where electron-phonon in- 
teractions lead to pairing. Fig. 3a is the analog of an inter- 
action mediated by a phonon propagating slowly relative 
to the fermions, at a velocity regulated by the Debye fre- 
quency, while Fig. 3b plays the role of a diagram using 
only the BCS model where the effective interaction is as- 
sumed to be instantaneous and attractive. Note that in 
principle this is a drastic approximation since the ratio 
of typical fermionic to phononic velocities is about 100 
in the low temperature superconductors. Nevertheless, 
in this context it is known that the BCS model captures 
the essence of the electron-phonon interaction, in partic- 
ular the s-wave symmetry of the superconducting con- 
densate. Thus, it is reasonable to believe that a similar 
situation will occur in models for the Cu-oxides based on 
the interaction between holes and spin-waves where the 
relevant ratio of velocities is likely of order 1 (since both 
are dominated by J in the low hole density limit, as ex- 
plained before). In other words, assuming that the inter- 
action mediated by the spin-excitations is instantaneous 
is an approximation that preserves the relevant quali- 
tative features of the problem, as the BCS model does 
in the electron-phonon problem. In this respect is that 
models of high-Tc such as those proposed in Ref. [0, 
where it is claimed that a peak in the density of states 
of strong correlation origin causes the appearance of an 
"optimal" doping, can be qualitatively correct and useful 
as a starting point to understand spin-wave mediated su- 
perconductivity. However, for a quantitative analysis the 
full spin-hole problem with retardation effects included 
must be taken into account, and that should be the main 
message extracted from the numerical studies discussed 
here. 



D. Implications for SO(5) scenarios 

Note that the results of this section have implica- 
tions for the recently proposed SO (5) scenario for the 
cuprates |^^. Once again if a snapshot of the system 
is taken at the intermediate time (dashed line) of Fig. 3a 
the holes will be found to be in a triplet state of mo- 
mentum q, and this state will have low-energy. However, 
as discussed before even though the lower energy spin- 
waves appear at q = Q, their number is maximized at 
q ~ (7''/2, 7r/2) and, thus, no particular value of the mo- 
mentum is a priori more important than others. This 
conclusion is in disagreement with the results of Ref. |2^] 
that favor q = Q over other momenta. In Fig. 4, the en- 
ergies of the two-hole ground state as a function of mo- 
menta for the cases of a total spin zero and one and using 
the t — J model are given. The results are exactly calcu- 
lated using a 18 site square cluster. In the triplet branch 
there is no major difference between the different mo- 
menta, in agreement with our discussion. The particular 
value q — Q is not specially important in a broad range 



of couplings from J /t = 0.2 to 1.6. The present reasoning 
is simple and avoids the use of complicated approximate 
symmetries of the t— J and Hubbard models pQ] . 




FIG. 4. Energy As=i(q) of the state with the lowest en- 
ergy in the subspace of two holes, momentum q, spin one, and 
parametric with J /t. The energies are referred to the overall 
ground state energy of two holes, which is a zero momentum 
spin singlet. The results are obtained exactly using an 18 site 
cluster. Results for 20 sites are very similar. 



III. EXACT DIAGONALIZATION 
A. Ladders 

Here it is addressed numerically whether the recent 
DMRG results for the behavior of two-holes in the 
ladder t — J model can be reproduced with other tech- 
niques. In particular, it would be desirable to obtain 
similar results in clusters accessible to Exact Diagonal- 
ization (ED) methods since the implementation of this 
algorithm is simple, calculations on a variety of i — J-like 
models can be carried out without much effort, periodic 
boundary conditions can be implemented, and dynamical 
properties can be analyzed. 

Using the ED method let us first investigate the rela- 
tive distance between holes in the case where two-hole 
bound states are expected. Fig. 5 shows exact results 
obtained on a 2 x 10 cluster with 2 holes and periodic 
boundary conditions. The probability P{d) of finding 
the holes at a relative distance d is presented against the 
coupling J /t. At every coupling, results are normalized 
such that J2d^('^) ~ ^- clear from the figure that 

for realistic values of J/t, such as 0.3 or 0.4, the highest 
chances indeed occur when d = -\/2 is the hole separation. 
Thus, DMRG and Lanczos techniques provide similar re- 
sults regarding this issue which is gratifying. However, 
Fig. 5 provides extra interesting information: P{^/2) is 
actually similar to P{1) and, thus, results in this context 
must be necessarily interpreted in a probabilistic sense 
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i.e. the distance with the highest chances is not neces- 
sarily much relevant for the problem. For instance, note 
that P{V^) ~ 0.28 which is substantially smaller than 1. 




FIG. 5. Exact Diagonalization results using a 2x 10 cluster 
with two holes. P{d) is the probability of finding the holes at a 
distance d apart (convention indicated in the figure) . Results 
for several couplings J/t are shown. 

This is better illustrated in Fig. 6 where P{d) at a fixed 
J/t is provided once one of the holes is fixed to an ar- 
bitrarily selected site of the 2-leg ladder. Working at 
J/t — 0.8 the second hole is mostly located on the chain 
the opposite to the first hole, more specifically in the 
three sites the closest to it. As J/t is reduced to 0.4, 
the second hole spreads further its wave function and the 
largest probability is at distance ^/2. However, visually 
it is clear that a better representation is to imagine the 
hole as moving freely within a small region in the vicinity 
of the fixed hole, rather than assign special importance 
to one of the possible distances in this region. In addi- 
tion, note that in the results of Fig. 5, and considering 
data for other couplings not shown explicitly, no abrupt 
changes were observed as a function of J/t for the hole 
distribution and smoothly the size of the hole pair wave 
function grows as J/t is reduced. Finally at J/t = 0.2 
the second hole seems delocalized on the whole ladder 
(suggesting the absence of binding, or a bound state of 
size larger than the cluster considered here). 

The DMRG studies isolating the configuration of 
holes where they are located along the diagonal of a pla- 
quette have also shown the existence of a strong spin 
singlet across the opposite diagonal of the same plaque- 
tte. A possible contribution to this result based on the 
d-character of the two-hole state was provided in Sec. II. 
As in the previous paragraphs, first it would be inter- 
esting to investigate if similar numerical information is 
obtained using ED techniques. In Fig. 7 results are pre- 
sented once again on a 2 x 10 cluster and with two 
holes. The bonds where the spin correlation (Si • Sj) has 
changed substantially compared with the undoped case 
are highlighted, following a convention similar as used in 
Rcf. |19|. At large J/t, and isolating the configuration 



where the two holes are in the same rung (since it has 
the highest chances) the spins along the nearest-neighbor 
rungs form a stronger spin singlet than in the undoped 
ladder. At J/t = 0.4, and now with two holes in the 
■\/2-configuration, the ED results confirm the previous 
DMRG calculations since indeed a robust spin singlet 
is found in the same plaquette where the holes are lo- 
cated. The considerable strength of the diagonal spin 
singlet naturally weakens the bonds that link the pla- 
quette being analyzed with the rest of the ladder. If J/t 
is further reduced to 0.2, now the most likely hole con- 
figuration have the holes at distance VS- Concentrating 
on such configuration, it is interesting to notice that here 
there are two spin singlets along the plaquette diagonals 
(see Fig. 7). Both are strong, causing the bonds between 
them to be weaker than in the undoped case. However, 
appealing as this picture may seem, the results in Fig. 6 
suggested that the holes are basically unbounded for this 
coupling. Thus, extracting conclusions out of snapshots 
of hole configurations is somewhat risky. 



J/t = 0.8 




J/t = 0.4 




J/t = 0.2 




FIG. 6. Ground state results obtained on a 2 x 10 cluster 
with two holes solved exactly. One hole is fixed at the position 
denoted by the open circle. The area of the gray circles is 
proportional to the probability of finding the other hole at a 
particular site. Results for several couplings are shown. 

As a partial summary, in this section it was found that 
(i) ED techniques reproduce the pattern of spin singlets 
found previously in DMRG studies once the position of 
holes are fixed, but (ii) several other hole configurations 
have comparable probabilities. The two hole bound state 
for realistic values of J/t has a finite size and it resembles 
a bi-spin-polaron, with two particles moving quasi-freely 
inside a region of space regulated by the coupling. There 
are no sharp hole distances dominating the ground state 
of the problem. The spins inside the two-hole wave func- 
tion form strong spin singlets which is natural since in 
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the undoped case the ground state has an RVB character. 
However, it was here argued that the d-wave properties 
of the state contribute at least in part to the strong spin 
singlet formation. 



J/t = 0.8 



J/t = 0.4 




FIG. 7. Results for two holes obtained at several couplings 
J/t on a 2 X 10 cluster solved exactly. The holes are fixed at 
the position with the highest chances in the ground state. 
Working in this subspace, bonds where the nearest-neighbor 
spin-spin correlation has the largest variation with respect to 
the undoped case are indicated. Solid (dashed) bonds indicate 
correlations which are larger (smaller) than in the undoped 
case by an amount larger than 20%. The thickness of the lines 
is proportional to the change observed. 

B. 2D clusters 

Calculations similar as those in the previous subsection 
are here repeated using finite two-dimensional clusters. 
Fig. 8 shows P{d) for a -\/20 x \/20 cluster solved exactly. 



0.6 




0.0 0.5 1.0 1.5 

J/t 



FIG. 8. Exact Diagonalization results using a 20 sites 
square cluster with two holes. P{d) is the probability of find- 
ing the holes at a distance d apart (convention indicated in 
the figure). Results for several couplings J/t are shown. 



As in the case of the 2-leg ladders, a finite window 
in J/t exists where the distance d = \f2 has the high- 
est chances, in agreement with Ref. However, if 
J/t = 0.3 or 0.4 is considered realistic, then holes at 
d — \/5 have the largest probability. However, as in the 
case of ladders there are several distances with similar 
weight i.e. P{d) is not sharply peaked at one particular 
value of d. Actually, Fig. 9 provides a pictorial represen- 
tation of the probability of having one hole at a given 
position once the other hole is fixed at an arbitrary site, 
at three values of J/t. As in the case of ladders, for large 
J/t the wave function of the second hole is localized in 
the immediate vicinity of the first one. Reducing J/t the 
extension of the wave function smoothly increases, and 
at J/t — 0.2 it covers the whole cluster suggesting the 
absence of a bound state. Thus, it seems that the bound 
state of two holes in 2-leg ladders and planes share many 
similarities, at least in the regime of small and interme- 
diate size pairs. This universality may be caused by the 
presence of robust antiferromagnetic correlations in both 
systems. The bound state main features seem indepen- 
dent of the long distance (gapped vs gapless) character- 
istics of the undoped ground state. 

The study of the pattern of spin singlets near holes 
shown in Sec.III.A can be repeated for the 2D clusters as 
well, and results are given in Fig. 10. At large J/t = 1.6, 
the holes are located with the highest chances at a dis- 
tance of one lattice spacing. As in the case of 2-leg lad- 
ders, the bonds parallel to the location of the holes con- 
tain strong spin singlets, which weakens the neighboring 
bonds along the same direction. When J/t is reduced 
eventually having the holes along a plaquette diagonal 
becomes the configuration with the highest chances (with 
the caveats of the previous paragraphs). In this case the 
spin singlet along the diagonal the opposite to where the 
holes are located is very strong, and as a consequence the 
other bonds associated with this spin singlet are weak, as 
it occurs for ladders. At small J/t = 0.2, once again two 
strong spin singlets are found between the holes which 
themselves have the highest chances of being located at 
distance -s/S, again similarly as observed for the same 
coupling on ladders. 

A possible explanation for the presence of strong sin- 
glets was provided in Sec. II where it was proposed that 
they arise as a natural consequence of the dx2-y2 symme- 
try of the two- hole ground state in a Neel background. To 
verify this statement, ED studies for the spin anisotropic 
t — J model were carried out. In this case the Heisenberg 
Si ■ Sj interaction is replaced by Sf + \{Sf + Sf Sf). 
If A = 1 (A = 0) the Heisenberg (Ising) limit is recovered. 
Working at J/t = 0.6 as an example, on 16 and 20 site 
clusters, and concentrating on two holes in the d^^-y^ 
subspace |2^ at distance \/2, the numerical results show 
that as A is reduced from 1 to the {SfS^) correlation 
(along the diagonal the opposite to where the holes are) 
is negative, of similar value in this interval of A (e.g. on 
the 20 site cluster, (S'f 5'p 0.21 at A = and 0.19 
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FIG. 9. Ground state results obtained on a 20 sites square 
cluster with two holes at J/t — 1.6, 0.4, and 0.2 solved exactly. 
One hole is fixed at the position denoted by the open circle. 
The area of the gray circles is proportional to the probability 
of finding the other hole at a particular site. More than 20 
sites in the clusters are shown for clarity (periodic boundary 
conditions were used). 

at A = 1). The value of the correlation changes smoothly 
as a function of A implying that at least in part the pres- 
ence of strong diagonal singlets in the Heisenberg limit is 
indeed caused by effects contained in Neel backgrounds, 
as explained in Sec. II. Very similar results have also been 
obtained for the t — J, model on a 2 x 8 ladder: consid- 



ering J/t = 0.4 as example, (5? 5?) 0.21 at A = 1.0 

J/t =1.6 



i8i 



J/t = 0.8 



mi 



J/t = 0.2 




FIG. 10. Results for two holes obtained at several cou- 
plings J/t on a 20 sites square cluster using the Exact Diago- 
nalization technique. The holes are fixed at the position with 
the highest chances in the ground state. Working in this sub- 
space, bonds where the nearest-neighbor spin-spin correlation 
has the largest variation with respect to the undoped case are 
indicated. Solid (dashed) bonds indicate correlations which 
are larger (smaller) than in the undoped case by an amount 
larger than 20%. The thickness of the lines is proportional to 
the change observed. More than 20 sites in the clusters are 
shown for clarity (periodic boundary conditions were used). 
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and ~ —0.22 at A = 0.4. Then, the effect appears also in 
ladders and in this case it complements the tendency to 
form singlets which is natural in this type of geometry. In 
addition, similar studies for the case A = have allowed 
us to verify that in the subspace of s-wave states the 
diagonal singlet is now a triplet, also in agreement with 
the predictions of Sec. II. Then, this establishes that the 
antiferromagnetic fluctuations and the d-wave character 
of the two-hole bound state contribute significantly to the 
formation of the diagonal strong spin singlet in ladders 
and planes. 



IV. TRUNCATED LANCZOS 

It would be interesting to check if the results obtained 
in the previous section survive an increase in the lattice 
size, specially for the case of the two-dimensional clus- 
ters. A priori the strong similarities between the DMRG 
and ED results on small clusters for the case of 2-leg lad- 
ders suggest that the two-hole bound state is not much 
affected by size effects, at least for the case of tight pairs. 
Nevertheless, even if only for completeness, the issue of 
size effects is here addressed using the so-called "Trun- 
cated" Lanczos (TL) approach [g2|. This technique is 
currently under development and it has the advantage of 
sharing the same good features of the ED method, spe- 
cially (i) the ability of producing dynamical information, 
(ii) the generation of results using states with momen- 
tum as a good quantum number (e.g. if periodic bound- 
ary conditions are used), and (iii) the ability to study 
Hamiltonians that include interactions at intermediate 
and large distances if needed. Thus, in addition to verify- 
ing the results previously discussed, this section has as an 
extra goal the test of the TL method in a physically rel- 
evant case. Although it is doubtful that the TL method 
will reach the same accuracy of the DMRG approach in 
the study of static properties of quasi-one-dimensional 
systems, it may provide an intermediate algorithm be- 
tween ED and DMRG keeping some of the main features 
of both techniques. Actually, recent developments |2^ in 
this context using an exact change of basis to work with 
better degrees of freedom for the truncation procedure 
have already produced promising results that may trans- 
form the TL method into a more widely used technique 
for the study of correlated electrons. 

Since the details of the Truncated Lanczos approach 
were described before in the literature p^ , here only a 
brief summary will be included. In the first step of this 
method one state of the full basis (belonging to the sub- 
space of momentum and total spin z-projection that will 
be investigated) is selected to start the iterations. In the 
analysis below the zero momentum Neel state was used. 
Note that this starting point is somewhat inefficient for 
the case of ladder systems where the ground state has no 
long-range order and a spin-gap, and the recent develop- 
ments 1 28 mentioned before have indicated that a better 



starting point would be to use rung spin singlets in the 
singlet-triplet basis. Nevertheless, for the particular case 
of two holes in the t — J model the approach in the 5^- 
basis seems accurate (as shown below) and, thus, results 
in the new basis will be postponed for a future publica- 
tion. After the initial state is chosen, the Hamiltonian 
is applied several times producing a basis set of a few 
hundred states. In this space, generated dynamically by 
the Hamiltonian, a Lanczos diagonalization is performed. 
The ground state wave function is analyzed and only the 
basis states with a weight |cp > A are kept. The cutoff A 
is selected such that the number of states after the trun- 
cation procedure is performed is only about 50% or less 
of the size of the matrix being studied. The procedure 
is repeated several times i.e. growing the space, diago- 
nalizing in the generated basis, truncating to a fraction 
of it ("back and forth" procedure), until the available 
amount of memory is exhausted. This slow-growth ap- 
proach has proven to work very well in some cases such 
as the t — Jz model |Q , and it is believed to produce rea- 
sonable results in spin-gapped models. However, below 
it will be shown that it works also for two-dimensional 
clusters that have low-energy excitations. 



A. Truncated Lanczos on Ladders 

Fig. 11 shows the ground state energy as a function of 
the size of the basis for a 2 x 16 cluster with 2 holes. This 
particular example can not be solved exactly with present 
day computers since the total Hilbert space contains ^ 
10^ states, even after translational invariance is used to 
reduce the basis size. However, Fig. 11 suggests that the 
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J/t = 0.4 

J/t = 0.2 
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-6.5 



1000000 2000000 
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FIG. 11. Ground state energy obtained using the Trun- 
cated Lanczos approach working at J/t = 0.2 and 0.4, as a 
function of the size of the basis used. The cluster has 2 x 16 
sites and there are two holes. The oscillations are produced 
by the "back and forth" procedure described in the text. For 
both couplings the energy is measured with respect to the en- 
ergy of the Neel state with two holes in the same rung, namely 
E = -J(48- 5)/2. 
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TL approach with only ~ 2 x 10^ states produces an 
energy converged up to 3 significant figures. The oscilla- 
tions in the energy shown in Fig. 11 as the dimension of 
the Hilbert space grows are a consequence of the "back 
and forth" procedure described before. 

In the absence of an exact value for the ground state 
energy of the cluster considered in Fig.ll, it is possible 
to judge the accuracy of the approximate state obtained 
using the TL method by evaluating correlation functions. 
In Fig. 12 a calculation similar to that reported in Fig. 7 
is presented, namely fixing J /t the configuration of holes 
with the highest chances is selected, and in that sub- 
space the spin correlations in the vicinity of the holes 
are shown. The links where the nearest-neighbor spin- 
spin correlation functions have the largest variation with 
respect to the undoped system (also calculated with the 
TL method) are highlighted. Qualitatively the results are 
very similar to those found using the 2 x 10 cluster, i.e. 
for J/t = 0.4 the \/2 hole configuration has the highest 
chances and the spin singlet along the opposite diagonal 
is strong, while for J/t = 0.2 the holes are at distance 
•^5 with two spin singlets formed in between. The ex- 
cellent agreement between the results obtained with the 
TL method compared with ED and DMRG techniques 
is somewhat surprising since in the TL method applied 
to the basis the spin-spin correlations away from the 
holes indicate correlations stronger than expected for a 
spin-gapped system. However, in spite of this fact the 
behavior of holes is not affected, reinforcing the notion 
that the physics of tight hole pairs is dominated by the 
presence of spin correlations at short distances indepen- 
dently of their long-range behavior. Similar conclusions 




J/t = 0.4 



J/t = 0.2 



FIG. 12. Results for two holes obtained at several cou- 
plings J/t on a 2 X 16 cluster using the Truncated Lanzcos 
algorithm keeping ~ 2 x 10® states. The holes are fixed at 
the position with the highest chances in the ground state. 
Working in this subspace, bonds where the nearest-neighbor 
spin-spin correlation has the largest variation with respect to 
the undoped case are indicated. Solid (dashed) bonds indicate 
correlations which are larger (smaller) than in the undoped 
case by an amount larger than 20%. The thickness of the lines 
is proportional to the change observed. 



Fig. 13 shows P{d) for the 2 x 16 cluster, and a couple of 
couplings. Several distances have comparable probabili- 
ties. In particular for J/t — 0.2 there is an approximate 
plateau in P{d) between d = 1 and 3. 




FIG. 13. Probability P{d) of finding the holes at a dis- 
tance d on a 2 X 16 cluster with two holes, calculated using 
the Truncated Lanczos algorithm. Results for two couplings 
J/t are presented. 



B. Truncated Lanczos on 2D Clusters 

Encouraged by these results, the TL approach was ap- 
plied to the 2D V32 X \/32 cluster with 2 holes. Although 
the basis needed to solve exactly this problem can be re- 
duced by a factor 4 compared with the basis used for 
the ladder (due to rotational invariance on square clus- 
ters), the problem is difficult to solve exactly |^ and 
in addition it has not been addressed with DMRG tech- 
niques. Fig. 14 contains TL results using ~ 2 x 10^ states 
isolating the hole configuration with the highest chances 
and studying the spin background in its vicinity. Once 
again, the results are almost identical to those found us- 
ing smaller clusters (Fig. 10). The dynamically generated 
strong plaquette-diagonal spin singlet clearly appears in 
the 2D cluster study, similarly as it occurs in 2-leg lad- 
ders. Fig. 15 shows that the same agreement with smaller 
cluster calculations occurs once the hole-hole correlation 
(n(O)n(r)), where n(r) is the hole number operator at 
site r, is calculated. For J/t — 0.4 the correlation is 
maximized at distance -\/2. However, remember that a 
given site has 4 neighbors at distance -\/2 but 8 at dis- 
tance -s/S, and thus the maximum probability P{d) is at 
d = -\/5 in the case of J/t = 0.2. Fig. 16 shows similar in- 
formation in a different representation. The results are in 
good qualitative agreement with those found in smaller 
clusters. 



can be reached analyzing other observables. For instance, 
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J/t = 0.4 



FIG. 14. Results for two holes obtained at J/t = 0.4 on a 
32 sites square cluster using the Truncated Lanczos algorithm 
keeping ~ 2 x lO'' states. The holes are fixed at the position 
with the highest chances in the ground state. Working in this 
subspace, bonds where the nearest-neighbor spin-spin corre- 
lation has the largest variation with respect to the undoped 
case are indicated. Solid (dashed) bonds indicate correlations 
which are larger (smaller) than in the undoped case by an 
amount larger than 20%. The thickness of the lines is pro- 
portional to the change observed. 



0.10 
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FIG. 15. Hole-hole correlation {n(0)n(r)) vs distance r 
obtained with the Truncated Lanczos method on a 32 sites 
square cluster and two holes keeping ~ 2 x 10^ states. The 
couplings are indicated. 
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J/t = 0.4 



FIG. 16. Ground state results obtained on a 32 sites 
square cluster with two holes at J/t — 0.4 using the Trun- 
cated Lanczos algorithm keeping ~ 2 x lO'^ states. One hole 
is fixed at the position denoted by the open circle. The area 
of the gray circles is proportional to the probability of finding 
the other hole at a particular site. 



V. INFLUENCE OF HOLE-HOLE 
ELECTROSTATIC REPULSION 

Results such as those contained in Figs. 6, 9 and 16 
clearly imply that a pair of holes form a bound state 
on 2-leg ladders and planes in the realistic regime of 
couplings close to J/t ~ 0.4, in agreement with a vast 
amount of previous literature At this coupling the 
"size" of the pair (defined as the mean distance between 
holes) is about a couple of lattice spacings. Since for 
J/t ~ 0.2, the bound state seems to have disappeared 
(according again to the same Figures), then the pair size 
must change very rapidly in the window of couplings be- 
tween 0.2 and 0.4, where a "critical" value J/t\c must 
exist leading to hole binding. 

Let us analyze in more detail a particular case such as 
J/t — 0.4 which is located in a window of couplings "a 
priori" presumed to be realistic. Here holes are mainly 
located at distances 1, ^/2, and V5 lattice spacings from 
each other. Such small distances naturally raise a couple 
of concerns: first, experimental results for YBCO have 
suggested that its coherence length ^ is about 15 A which 
roughly corresponds to 4 lattice spacings ^ . Doped La- 
cuprates have an even higher ^ of about More recent 
estimations arrive to ^ = 18. 3A and 22. 7A for 1% Zn- 
doped underdoped YBCO and La-214, respectively [^ . 
These lengths are certainly much shorter than those ob- 
served in low temperature superconductors, but still do 
not locate the cuprates in the regime where ^ is smaller 
than the mean distance between carriers. If ^ is inter- 
preted, as usually done, as the size of the Cooper pairs 
in the superconducting condensate and if the t — J model 
is used, then J/t must be fine-tuned closer to its critical 
value J/t\c to have hole pairs of size equal to the exper- 
imentally measured ^. In other words, at J/t = 0.4 the 
hole pairs in the t — J model are substantially smaller in 
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size than needed to represent the cuprates. Then, from 
this point of view the pecuhar properties of a very tight 
bound state are not much relevant, and only the qualita- 
tive aspects of the problem can be extracted from studies 
of the pure t— J model at J /t — 0.4. 

In addition, having pairs of small size raises concerns 
regarding the stability of such pairs once the Coulomb 
interaction is considered. While the on-site electronic re- 
pulsion is usually taken into account in the framework of 
strongly correlated electrons, the repulsion at the next 
relevant distance of one lattice spacing is usually ne- 
glected (with an exception being the analysis of stripe for- 
mation in the cuprates pl|]). Other authors have raised 
similar concerns . A rapid estimation shows that the 
Coulombic electrostatic effect cannot be neglected since 
the bare potential energy between two charges at one lat- 
tice spacing from each other is Vnn = e^/a = 3.8 eV , 
where a = 3.8A as in the Cu-oxides compounds. Cer- 
tainly this repulsion is considerably larger than the at- 
traction expected between holes which is regulated by 
J ^ Q.leV . This source of attraction would be destroyed 
by the bare Coulomb interaction, unless retardation ef- 
fects are important. 

However, it is reasonable to consider a more optimistic 
scenario where the Coulomb repulsion at distance of one 
lattice is influenced by other orbitals in the Cu-ions, po- 
larization of oxygen [ |33| , and possible metallic screening 
in the doped regime. To gain some intuitive insight about 
the influence of these effects let us first simply divide 
Vnn by the dielectric constant Coo obtained experimen- 
tally for the high-Tc compounds (certainly being aware 
that the result will be qualitative at best since one is 
not supposed to use such a constant for the short dis- 
tance effects discussed in this section). Although there 
is a large range of estimations in the literature, a num- 
ber near ^ 30 — 40 is a reasonable assumption . 
However, even considering this factor the new estimation 
of the Coulomb repulsion still gives a number competing 
with the order of magnitude of the attraction caused by 
antiferromagnetism, namely Vnn = e^/(eooa) = 0.1 eV . 
A variety of other calculations lead to qualitatively simi- 
lar results. The Cu-0 repulsion (V^d) has been estimated 
to be ~ l.QeV [Q, i.e. smaller than the bare Coulomb 
interaction at distance 1.9A roughly by a factor 8. If this 
same factor is applied at distance of one lattice spacing 
the Cu-Cu repulsion {Vdd) should be about 0.5 eV. In ad- 
dition, note that the binding energy of a hole pair in two- 
dimensional clusters for J /t — 0.4 is only As ^ 0.2 J ^] 
i.e. while the hole attraction is regulated by J the numer- 
ical coefhcient that takes into account the fact that holes 
actually lose kinetic energy by forming the pair magnifies 
even more the damaging effects of the Coulombic repul- 
sion. 

Then, electrostatic effects will likely make unstable the 
very tight pairs observed in the computer simulations de- 
scribed before. To analyze explicitly this effect in Fig. 17 
results for the probability of finding the holes at a given 
distance d are provided for the case of a 2-leg ladder . 



0.3 




FIG. 17. Probability P{d) of finding the holes at a dis- 
tance d on a 2 X 10 cluster with two holes at J /t = 0.4, 
calculated using the Lanczos algorithm. Results for several 
values of the coupling V regulating a nearest-neighbor repul- 
sion between holes are presented. 

The Coulomb repulsion is included adding to the t — J 
Hamiltonian a term 1^ ^i'^j ^ where rii is the number 
operator at site i and the rest of the notation is stan- 
dard. The results show that the relevance of very small 
distances is actually lost already when V ^ 2J ^ 0.2eV. 
At larger values of the repulsion such as y = 4J the 
probability has a broad plateau between distances 2 and 
3, and the bound state seem to disappear when V is in- 
creased slightly further |^^. Note, however, that if the 
hole-hole instantaneous Coulombic repulsion is assumed 
to be short-ranged more extended bound states may sur- 
vive the inclusion of such electrostatic energies since their 
size increases as J/t is reduced to its critical value. But 
this once again will imply that a typical bound state size 
will be larger than just a couple of lattice spacings. 

There are two possible ways to avoid the problem of 
the Coulomb interaction at short distances. One is by 
considering retardation effects, and the other by finding 
other sources of screening that could reduce drastically 
the previous estimations of Vnn- Let us consider retar- 
dation first: the discussion in the previous sections have 
indeed suggested that this effect is important in the t — J 
model and it may avoid the problem of the Coulombic 
interaction similarly as it occurs in electron-phonon sys- 
tems, namely a hole could scramble the spin order in a 
region of space during its movement, and a second hole 
could take advantage of this distortion. If one hole fol- 
lows the other at a distance of a few lattice spacings, then 
the short-distance Coulomb repulsion can be avoided and 
the results would be compatible with estimations of ^ in 
the literature. As a simple illustration in Fig. 18, a hole 
creates a string excitation |3^ and a second hole heals 
the damage in the spin background at a later time [ p9| . 
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but keeping some distance to avoid repulsions. A similar 
process can potentially take place in the case of ladders 
as also shown schematically in Fig. 18. Considering the 
ground state as dominated by rung singlets, the move- 
ment of a hole produces the transformation of one of 
those singlets into a singlet along the diagonal of a pla- 
quette. Then, moving a hole introduces a damage in the 
spin configuration similar to what occurs with strings in 
the planes. The second hole can heal this damage as 
shown in Fig. 18. 



(a) 




(b) 



FIG. 18. (a) Possible behavior of holes in an antiferro- 
magnetic background that avoids the short distance Coulom- 
bic repulsion. The first hole moving, e.g., to the left creates a 
string of spins (highlighted in the figure) incorrectly aligned 
with respect to the background. The second hole can heal 
this damage following the first hole at some distance. The 
picture is possible if the spin excitations have a large lifetime 
such that the string of spins is not erased by spin fluctuations 
before it is used by the second hole to improve its energy; (b) 
Similar ideas as in (a) but now in the context of ladders. Here 
the spin background is not antiferromagnetic but the spins 
mainly form singlets along the rungs. A hole moving along 
the chains creates a string of diagonal singlets. The second 
hole can heal this damage, keeping some distance from the 
first to avoid the Coulombic repulsion. 

The second possible solution invokes a source of screen- 
ing not considered in the previous discussion. In the ab- 
sence of other holes and using the t — J model where 
electron-hole pairs cannot be created, metallic screen- 
ing effects should not affect the results and, thus, the 
Coulombic interaction will actually be of long-range. 
However, in the presence of a finite density of holes, 
a simple Thomas-Fermi approximation can be used to 
make a rough estimation of a possible metallic screen- 
ing length A. Considering a hole density n '-^ 7 y. 
IQ^^cm^^ and using standard textbooks equations 
for three dimensional metals, A is found to be ~ O.SA, 
which is small indeed and could potentially drastically 



reduce the effects of the Cu-Cu Coulomb interaction. 
However, it is clear that the Cu-oxides are very different 
from the three dimensional metals where Thomas-Fermi 
approximations are qualitatively reliable. Thus, the is- 
sue of whether metallic screening can reduce the effects 
of Coulomb repulsions between holes in the cuprates re- 
mains an open question that deserves further studies. 



VI. CONCLUSIONS 

In this paper the problem of two holes in a spin back- 
ground with robust antiferromagnetic correlations was 
revisited. Using Exact Diagonalization and Truncated 
Lanczos techniques applied to planes and ladders with 
up to 32 sites the typical distances between holes was 
studied. It was observed that the maximum probabil- 
ity occurs when the holes are at lattice spacings 
of each other, in agreement with previous calculations. 
However, other distances were found to be equally rel- 
evant. An intuitive explanation for these results in a 
real-space picture, as well as diagrammatically was pro- 
vided. Strong spin singlets are found near two holes in 
a d-wave state. An explanation for this effect based 
on a Neel background was presented, complementing 
other approaches ||l^ based on spin disordered back- 
grounds. It is concluded that holes in their movement 
create spin excitations with a non-negligible lifetime, 
and thus retardation effects are important in the t ~ J 
model. However, the instantaneous approximation cap- 
tures properly the qualitative aspects, specially the sym- 
metry of the bound states which is in the dx^-y^ chan- 
nel. The short size of the pairs raises concerns regarding 
the stability of the bound states if NN Coulomb inter- 
actions are included. Estimations of the strength Vnn 
of these interactions were here provided. A numerical 
study showed that the bound states increase their size 
as Vnn grows, and the typical hole distance in the pair 
can become comparable to the experimentally measured 
coherence length, which is about 4 lattice spacings, even 
if the values of J /t used produced unrealistic small pairs 
in the absence of the NN repulsion. 
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